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In this work we present a molecular dynamics simulation of a FFM experiment. The tip-sample 
interaction is studied by varying the normal force in the tip and the temperature of the surface. The 
friction force, cA, at zero load and the friction coefficient, fj,, were obtained. Our results strongly 
support the idea that the effective contact area, A, decreases with increasing temperature and the 
friction coefficient presents a clear signature of the premelting process of the surface. 
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I. INTRODUCTION 



Friction is one of the oldest phenomenon studied in 
natural sciences. In a macroscopic scale it is known 
that the friction force between surfaces satisfies the fol- 
lowing rules: (1) The friction is independent of contact 
area between surfaces; (2) It is proportional to the nor- 
mal force applied between surfaces and (3) The kinetic 
friction force is independent of relative speed between 
surfaces lj. Considering that friction is the result of many 
microscopic interactions between the building atoms at 
the surfaces, it must depend on factors as roughness , 
temperature and the energy dissipation mechanism at 
the surfaces. Therefore, to understand its macroscopic 
behavior it is necessary to understand in details the dy- 
namics of interaction between atoms in the surfaces in 
contact. In 1987, C. M. Mat et al 2] have used, for the 
first time, the Friction Force Microscope (FFM) to in- 
vestigate friction in nano-scale. That kind of microscope 
allows the experimentalist to produce essentially a single 
contact between a sharp tip, of atomic dimensions, and 
the surface. U I3- Its basic principle is shown in fig- 




FIG. 1: Schematic view of a FFM apparatus. 151 

ure^ In this apparatus the tip stays in touch with the 



surface while it moves at constant velocity, v, parallel to 
the surface. The resulting force in the cantilever makes 
a torsion on it. This torsion can be observed by optical 
techniques. One of the more striking effects observed is 
the stick-slip phenomenon: The friction force plotted as 
a function of time exhibits a sawtooth behavior. 
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Analytically, the motion of the tip can be modelled as 
follows. Forces in the tip are a sum of two terms: An 
interaction force between the tip and the surface due 
to the periodic atomic arrangement of the lattice and 
a force due to the cantilever. Some authors by using 
this approach were able to reproduce several features 
of the friction at a nanoscopic scale [a, @ ■ In this work 
we use molecular dynamics (MD) simulation to study 
the friction phenomenon at the atomic scale. In figure 
121 we show a schematic view of the model we have used 
in our simulation to reproduce the FFM mechanism 
(Figure ^| . The tip is represented by a single particle 
that interacts with the surface through a convenient 
potential. The springs represent the mechanism we 
have used to vary the normal force (z direction) and 
to measure the lateral force (x and y directions). By 
measuring both forces it is possible to study the friction 
force behavior under several circumstances. 



In a recent work Resende and Costa [lfj using molecu- 
lar dynamic simulation have studied the migration of an 
individual atom on the surface of a 12 — 6 Lennard- Jones 
bec crystal. They argued that an observed anomaly oc- 
curring in the diffusion constant can be the signature of a 
pre-melting process. The migration of an ad-atom at the 
surface may occurs by three mechanisms. At low tem- 
perature the adsorbed particle can move through chan- 
nels on the surface since thermal motion of atoms at the 
surface have low amplitude. Once temperature rises we 
reach an intermediate state. The surface starts to melt 
so that the channels are closed and ad-atoms are stuck in 
the vicinity of a surface atom. The situation persists until 
the ad-atom is thermally activated and random- walk dif- 
fusion occurs. In summary, the diffusion constant should 
present a minimum at the intermediate region. Under the 
point of view of friction we may ask what is the effect of 
this phenomenon over friction. For two macroscopic slid- 
ing surfaces we may not expect to distinguish the first 
two process since the contact area is large compared to 



interatomic distance. As temperature rises the surface 
is lubricated by melted atoms, we may expect a smaller 
friction coefficient. The situation is quite different for a 
small tip in contact with the surface. 




FIG. 2: Schematic view of our Friction Force apparatus. 

In the following we describe a series of molecular dy- 
namics computer simulation of the interaction of a small 
tip with a surface. Computer simulations give us a 
very convenient instrument to study tribological pro- 
cesses. They allow controlled "numerical experiments" 
where the geometry, sliding conditions and interactions 
between atoms can be varied at will to explore their ef- 
fect on friction, lubrication, wear and to reach unacces- 
sible regions for experimentalists. Unlike laboratory ex- 
periments, computer simulations enable us to follow and 
analyze the full dynamics of all atoms. A wide range 
of potentials have been employed in numerical studies of 
tribology. For example, simulations of metals frequently 
use the embedded atom method, while studies of hydro- 
carbons use potentials that include bond-stretching and 
so on. We will not concentrate ourselves in any specific 
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material rather in general aspects of the problem under 
consideration. Due to the general character of our study 
we content ourselves by using the Lennard- Jones (6 — 12) 
potential which is commonly used in studies of general 
behavior. We will consider that the interaction of the 
tip with the surface does not disturb very much the elec- 
tronic structure of the system. This consideration means 
that we do not have to deal with the quantum aspects 
of the problem. This simplification let us save a lot of 
computer time. If, for one side, we lose details of the 
considered phenomenon, we gain in that we can study 
true dynamical and temperature dependent models. 

This work is organized as follows. In the section ITT1 we 
introduce general aspects of the numerical method used, 
in the section ITTT1 we present our results and in section llVl 
we discuss and present some conclusions. 



II. SIMULATION BACKGROUND 

Our simulation is carried out by using molecular dy- 
namics (MD) simulation. A schematic view of the sim- 
ulation arrangement we have used is shown in figure |21 
Three springs of elastic constants k x , k y and k z are at- 
tached to the particle M that represents a tip point. This 
arrangement allow us to measure normal (F z ) and paral- 
lel (F x , F y ) forces on M. The surface is represented by an 
arrangement of particles which interact with each other 
and with the mass M through a truncated Lennard- Jones 
(6 - 12) (LJ) potential 
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where 4>i,j{ri.j) is the complete LJ potential, 
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The indexes i and j stands for position vectors rj and 
rj respectively, and 1 < i < N, where N is the to- 
tal number of particles and rij — \rj — J"j|. A cutoff, 
r c , is introduced in the potential in order to acceler- 
ate the simulation. If the force on a particle is found 
by summing contributions from all particles acting upon 
it, then this truncation limits the computation time to 
an amount proportional to the total number of particles 
N. Of course, this truncation introduces discontinuities 
both in the potential and the force. To smooth these dis- 
continuities we introduce the constant term (j>(r c ). An- 
other term (d(f>ij(ri t j)/dri t j) is introduced to re- 



move the force discontinuity. Particles in the simulation 
move according to Newton's law of motion, which gener- 
ates a set of 3A^ coupled equations of motion which are 
solved by increasing forward in time the physical state 
of the system in small time steps of size St. The result- 
ing equations are solved by using Beeman's method of 
integration[ll], [I3, ll3L Il4| . The system is arranged in 
4 layers with free boundary conditions in all directions. 
The first layer is frozen in a regular arrangement as in 
the (001) surface of a Lennard- Jones bec crystal in order 
to maintain the whole structure as flat as possible. 

With the tip far away from the surface we thermal- 
ize the system at temperature T. After thcrmalization, 
the tip is pushed in a direction parallel to the surface at 
constant velocity v p . For each simulation the distance 
between the spring and the surface is fixed at the start, 
so that we can control the perpendicular force on the tip. 
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FIG. 3: Energy as a function of temperature. The melting tempera- 
ture is estimated as the inflexion point, being around T m S3 l.le/k^ 



By measuring the size variation of the springs we can 
calculate the lateral, F x , and the perpendicular force, 
F z , on the tip. The temperature, T, of the surface can 
be controlled by using a velocity renormalization scheme 
(See for example |l5| and references therein). From the 
equipartition theorem we can write that 



( V 2 ) = 3 ^B T 
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(3) 



By controlling the value of (v 2 ) we can reach a chosen 
temperature Tf. An appropriated way to do that is by 
successive approximations. We chose a set of velocities 
{v}o so that we get (v 2 )o- By multiplying each velocity 
by a factor ap defined as 
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a first approximation to Tf is done. By evolving in time 
the system we can create sequences, T n , {v} n and {a} n , 
such that after a finite number of time steps the temper- 
ature of the system converges to T n sa Tf. The friction 
coefficient is calculated as the quotient 
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dF x 
dF z 



(5) 



Before we start the simulation we have to have an es- 
timative of the melting temperature, T TO , of the system. 
This is done by performing a preliminary simulation of 
the substrate. In Figure |3 we show the total energy per 
particle, E, as a function of temperature. The melting 
temperature is estimated as the inflection point of the 
curve. We find TL w 1.1 in accordance with earlier calcu- 
lations [nj, Ha, Ha, D3 • The velocity, position and forces 
are stored at each time step for further analysis. We 
measure the time t, temperature T and forces in units of 
ay^m/e, e/ks and e/a respectively. 



III. RESULTS 

We have simulated the FFM system for several tem- 
peratures and initial distances of the tip to the substrate 
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FIG. 4: The friction force, (F x ) as a function of normal force, (F z ) 
is shown for several temperatures. The (F x ) and (F z ) forces are 
measured in units of e/a. The figures, from a to e are for several 
different values of T = 0.25,0.44,067,085,1.05 respectively. The 
circles are the MD results and the straight line correspond to an 
adjust. 



or equivalently, the normal force in the tip. In figure 01 
we show a plot of our MD simulation results for the fric- 
tion force as a function of normal force for several tem- 
peratures. The Amonton's Law of friction states that 
frictional forces are proportional to the normal force and 
independent of the contact area. This type of behav- 
ior was observed in some systems by many authors, who 
fitted {F x ) to a linear function of both load,(F 2 ), and 
contact area, A: 



(F x ) = fi{F z ) + cA. 



(6) 



Here \i is the friction coefficient and the second term 
cA is interpreted as the friction force for zero normal 
force. Ringlein et. al.[lq showed that the friction forces 
in nanometric adhesive systems depends on the contact 
area, breaking down the Amonton's law. Because we use 
the LJ force between the tip and the surface is adhesive. 
In the following we present our results that strongly sug- 
gest that the Amonton's laws [lei 119] is violated when the 
friction force is considered as a function of temperature. 
In the figure OJleft) we show a plot of cA as a function 
of temperature. We can see that when the temperature 
increases the contact area or adhesion forces decrease. 
This behavior can be related to the fact that at low tem- 
perature the atoms at the surface perform low amplitude 
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FIG. 5: Plot of cA (left) and /i (right) as a function of Tempera- 
ture. The line is only a guide to the eyes. 
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FIG. 7: Histogram for the residence time. The normal force is for 
F z = -1.09, -1.02, -0.95, -0.42, -0.14. Beam size is tbeam = 5. 



-10 10 



FIG. 6: Path of the tip over the surface (XY plane) for sev- 
eral temperatures and normal forces. From a) to e) we have 
T = 0.25, 0.44, 067, 085, 1.05 respectively. The normal forces arc 
defined in fig |1] The plots are shown dislocated by a constant 
value in the y direction as a matter of clarity. 



jumps so that the number of collisions with the tip is low. 
In this case the effective contact area is high because the 
tip stays a long time close to the surface. However, when 
the temperature grows the number of high energy fluctu- 
ations of particles at the surface increases with a conse- 
quent increase in the number of collisions of high energy 
with the tip, decreasing the effective contact area. We 
also observe that the friction coefficient (Shown in figure 
Upright)) grows abruptly at T ~ 0.7e/fcfc. In the follow- 
ing we discuss the fact that this behavior can be related 
to the pre-melting of the surface. We show in fig. H3 a 
plot of the path of the tip over the surface for several 
temperatures and normal forces. As should be expected 
the paths are well defined for low temperatures becoming 
random as temperature grows. To understand what is 
going on, we have calculated the residence time of the tip 
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FIG. 8: Histogram for the residence time. The normal force is for 
F z = 2.05, 2.09, 2.13, 2.52, 2.56. Beam size is t beam = 5. 



in each site, defined as the time that the tip spent in the 
neighborhood of a specific site, i.e., its distance to some 
particular site being smaller than some reference distance 
S. With no loss of generality we choose 8 — a, the lattice 
parameter. In figure |7| and [3] we show the histograms 
for normal forces F z = -1.09, -1.02, -0.95, -0.42, -0.14 
and F z = 2.05,2.09,2.13,2.52,2.56 respectively for sev- 
eral temperatures. For negative values of F z the residence 
time is well defined even at high temperatures having its 
average at St « 50. The tip is immersed in the surface, 
so that, it can easily travel along channels on the surface 
of the crystal. In the figure |S1 we show the histograms 
for positive values of F z . The temperatures are the same 
as in figure ^ At low T they show a similar behavior 
as that for F z < 0. However, at higher temperatures, 
the residence time spreads out to the higher t region. 
We interpret this as a closing of the channels discussed 
above due to thermic motion of the particles at the sur- 
face. When temperature increases the particle gets more 
energy, which is eventually enough to push it from any 
specific neighborhood. 



IV. CONCLUSION 

We have performed a molecular dynamics simulation of 
a FFM experiment. Our results were obtained by vary- 
ing the normal force in the tip and the temperature of the 
surface. The behavior of the cA term in the Amonton's 
law (eq. 

EJ) and the friction coefficient were found to depend 
on the temperature. The cA term which measures es- 
sentially the effective contact area, A, between the tip 
and the surface were found to decrease with increasing 
T. The friction coefficient presents a sudden jump which 
seems to be related to the premelting processes of the 
surface. 
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